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ABSTRACT 

The  longitudinal  differential  equation  of 
motion  has  been  used  to  investigate  various 
aerodvnamic  expansion  techniques.  The  total  drag 
coefficient  was  expanded  using  conventional 
polynomials  and  splines  with  and  without  floating 
knot  locations.  This  paper  discusses  the  various 
techniques  and  approaches,  compares  results 
obtained  from  s imul taneous 1 v  fitting  four 
separate  flights  (time  vs  distance  measurements) 
and  outlines  the  potential  advantages  and/or 
disadvantages  of  the  various  aerodvnamic 
expansion  techniques.  It  is  believed  that  this 
is  the  fir-c  time  splines  have  been  used  in  the 
aerodvnamic  coefficient  estimation  process  and 
chat  these  results  and  techniques  are  germane  to 
other  applications. 

NOMENCLATURE 

A  Reference  area 

a  Coefficient  in  Equacion  1 

C;  Elopes  of  spline  segments 

'see  Equations  4,  S,  and  h) 

c p  T''tal  drag  coefficient 

l"r>n  Eero  angle  of  attack  drag  coefficient 

~U->  Second  order  drag  term 

'sop  Equations  ?  and  1' 

1  D',  Fourth  order  drag  term  (see  Equation  1) 

1  %  brag  variation  due  to  velocitv  change 

'I  Exponential 

m  Model  mass 

Number  of  straight  line  segments 
1  see  Equation 

'n  of  residuals  squared 


I  lie  |i.i|i.  r  i.  ill  .  ljri-,1  j  ,|,t  (  % 

<  ■ 1 1  * i  r M n"  nt  .in 1 1  ihiTrlnrr  I'  in  (hi-  puhlu  dunum. 


Si  Switches  for  spline  segments  \ 

(see  Equations  4.  ■>,  and  El  ' 

V  Velocitv  along  the  X  axis 

Vi  Instantaneous  velocitv 

Vref  Reference  velocitv 

X  Dovnrange  axis 

‘i  Instantaneous  pitch  angle 

i  Instantaneous  vaw  angle 

Total  instantaneous  angle  of  attack 

'i  Knot  locations  in  Equations  4,  5,  and  E 

o  Air  densitv 


Superscripts 
•  Time  derivative 

INTRODUCTION 

Prior  to  19E9,  the  prevalent  method  of 
analyzing  ballistic  spark  range  data  was  based  on 
the  linear  approximation  method  known  as  "linear 
theory"  developed  bv  Murphv*~^  and  others7 
Stated  briefly,  the  method  uses  a  closed-form 
solution  to  the  differential  equations  of 
motion.  This  approximate  solution  results  from 
assuming  a  linearized  aerodvnamic  model  where  t"' 
aerodynamic  force  and  moment  derivatives  are 
constant  with  angle  of  attack  (hence  the  name 
"linear  theory").  Murphv  extended  this  technique 
to  include  a  quasi-nonl  inear  analysis*’  where 
the  linear  aerodynamic  force  and  moment 
derivatives  are  reduced  in  such  a  manner  that 
certain  nonlinearities  could  he  obtained.  This 
quas i-nonl inear  analysis  requires  an  assumed 
functional  form  of  the  nonlinearitv  (normally  * 
quadratic  or  cubic  polynomial). 

In  19E9  Chapman  and  Kirk,  in  analyzing  fri 
'flight  data7,  documented  the  application  or  a 
technique  thev  called  parametric  d  i  f ferent i at > 
which  permitted  the  free  flight  differential 
equation  of  motion  to  he  used  directlv  in  t  r 
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a  correlation  process.  This  technique 
liminated  the  requirement  for  closed-form 
*  lotions  to  the  equations  of  motion.  However, 

.  js  still  required  to  assume  a  form  of  the 
nonl inear  it ies  in  the  equations  of  motion, 
general lv  these  forms  have  also  been  assumed  to 
h|1  polynomial  expansions  of  the  aerodynamic  force 
3nd  moment  derivatives  with  angle  of  attack8_q. 

During  the  past  several  vears,  data  analysts 
nave  discussed  the  possibility  of  using 
mathematical  splines  (two  or  more  mathematical 
^xr,rPS9 i ons  attached  end  to  end!  for  the 
coefficient  expansions.  These  splines  permit  the 
Slopes  of  the  aerodynamic  expansions  to  be 
discontinuous  and  would  offer  the  analyst  a  more 
reneral  aerodynamic  mode'.,  thereby  relieving  some 
of  the  requirements  of  assuming  the  form  of  the 
nonl inear ic ies .  This  paper  discusses  various 
coefficient  expansion  techniques  and  compares 
results  obtained  using  the  various  expansions. 

METHOD  OF  APPROACH 

In  order  to  evaluate  the  various  expansion 
techniques  (continuous  function  vs  splines),  we 
will  restrict  our  attention  to  a  simple  single 
degree  of  freedom  system  rather  than  the  full 
six-degree  of  freedom  svstem  described  in  various 
references. 


This  paper  will  examine  the  determination  of 
the  total  drag  coefficient  (Cp)  as  a  function 
of  instantaneous  angle  of  attack  and  velocity 
from  the  longitudinal  momentum  equation  and  the 
associated  experimental  measurements  of  distance 
traveled  vs  time.  The  differential  equation 
governing  the  longitudinal  momentum  is 

•  _  "  P  •  - 

tr 


“here  -  is  the  air  densitv,  A  is  the  bodv 
reference  area,  m  is  the  projectile  mass,  X  is 
the  longitudinal  down  range  distance,  and  (•) 
and  (••)  indicate  the  first  and  second 
derivatives  with  respect  to  time.  The  total  drag 
coefficient,  Cp  (,V),  depends  on  the  instan- 
taneous  total  angle  of  attack,  “Jsin-  .♦  sin-£\ 
and  velocity,  V  «  X.  Equation  'l)  assumes  that 
the  angle  between  the  velocity  vector  and  the  X 
axis  is  small . 


fxpansions  for  the  Total  Drag  Coefficient 


Several  expansion  techniques  for  Cp  were 
investigated.  The  first  two  involved  continuous 
functions,  beginning  with  the  classical  quadratic 
t'«r*ndence  on  angle  of  attack,  or 


filiation  (2)  represents  the  classical  expression 
or  which  MurphvWZ  developed  a  methodology  of 
'"draining  Cp^  and  Cp^  bv  plotting  the 
*ff*ctive  measured  drag  coefficient!  va  the  mean 


of  «2.  A  straight  line  through  this  data 
yields  the  intercept,  Cp  ,  and  the  slope, 

Cpj,  The  second  continuous  function  is  much 
more  versatile  and  is  valid  for  a  wider-  ranee  o' 
angle  of  attack  dependence  and  a  linear  vi  ;  >.•  in¬ 
dependence,  that  is 

i  ..  ‘  .  p  •  •  '  '  > 

Note  this  expansion  has  five  unknown 
coefficients,  Cp^.  a,  Cp,,,  Cp^  ,  an! 

Cpv  as  do  all  the  remaining  expansion 
techniques  considered  within  this  pappr.  The 
d3-  term  is  somewhat  unconventional  but  allows 
a  nonzero  slope  at  zero  angle  of  attack.  The 
Cp  coefficient  is  normally  small  but  ror 
high  drag  configurations  which  experience  large 
velocity  decavs  during  the  flight  or  when 
time-distance  data  obtained  from  several  flights 
(slightly  different  launch  velocities)  are 
simultaneously  analyzed  this  term  can  be  impor¬ 
tant.  This  Cp _  term  accounts  for  variations 
in  drag  coefficients  with  Mach  number  and 
Revnolds  number.  The  two  effects  cannot  be  sim¬ 
ply  separated  because  they  both  depend  linearlv 
on  velocity.  This  term  appears  in  all  -f  the 
expansions  discussed  except  for  the  clas'-ical 
quadratic  dependence  shown  in  Equation  '2K 

Several  expansions  using  spline  functions 
were  also  evaluated.  The  first  of  these  uses 
multiple  straight  line  segments  (see  Figure  1) 
and  the  general  expansion  is 


+  c.  < :  .  ,  ■*  t  i '  - 

i  i* !  :  1' 

where  i  is  the  index  for  the  spgoent.  '  s  ar.- 
the  knot  locations.  C;'s  are  the  slope-  or  oac! 
segment,  and  S;'s  are  the  switches  ' Sj  =  n  or 
l  depending  on  whether  or  not  the  instantaneous 
value  of  ■  is  within  the  range  of  the 
segment).  For  the  investigation  discussed,  the 
evaluations  of  the  various  expansion  techniques 
were  restricted  to  four  unknown  coefficients  plus 
the  Cp  term.  Hence  onlv  two  cases  utilizing 
straight  line  segments  were  considered.  The  first 
case  is  three  segments  (n*3)  with  the  knots  fixed 
at  j,  and  p.  These  knot  locations  were 
chosen  bv  dividing  the  range  into  three  equal 
parts,  V  max'3.  : 2*2  max max' 

The  five  unknown  coefficients  then  become 
CD0»  cl*  c2’  c3*  and  CD  •  Here  it 
should  be  noted  that  other  methods  of  dividing 
the  S  space  into  segments  and  selecting  the  knot 
locations  were  considered.  Initially,  it  was 
felt  that  the  segments  should  be  chosen  such  that 
Chose  associated  with  the  higher  angles  of  attack 
would  be  small  compared  to  the  segments  associ¬ 
ated  with  Che  smaller  angles  of  attack.  The 
reaaoning  for  this  was  that  it  was  assumed  that 
Che  rate  of  change  of  the  drag  coefficient  with 
respect  to  '  wss  much  higher  at  large  angles  of 
attack  thereby  requiring  smaller  segments. 
Although  this  assumption  is  certainly  true  for 
most  free  flight  range  data,  the  nature  of  a  well 
behaved  dynamically  stable  configuration  in  free 
flight  is  that  the  larger  initial  angles  of 
attack  rapidly  decrease  (damps)  during  the 
flight.  Hence  onlv  a  relatively  few  data  points 
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representing  the  initial  Mgh  angles  of  attack 
are  normal ly  obtained  compared  with  the  number  of 
data  points  associated  with  the  smaller  angles  of 
attack  (for  example,  see  Figure  2).  Considering 
this,  if  a  small  segment  was  chosen  for  the 
higher  angles  of  attack,  only  a  very  few  data 
points  would  fall  within  this  segment  thereby 
invalidating  the  resultant  slope  parameter. 

After  this  anomaly  was  recognized  it  was  felt 
that  perhaps  Che  knot  locations  should  be 
selecred  such  that  each  segment  contained  an 
equal  number  of  data  points.  However  it 
insediately  became  obvious  that  due  to  the  nature 
of  the  data,  a  very  small  segment  resulted  at  the 
smaller  angles  of  attack,  where  it  wasn't  needed, 
and  a  large  segment  resulted  at  the  large  angles 
of  attack  where  a  small  segment  was  desired  (see 
Figure  1).  With  these  considerations  is  mind, 
dividing  the  •'  range  into  equal  parts  appeared  to 
be  a  reasonable  compromise.  However,  if  this 
technique  is  applied  to  other  free  flight  data 
fi.e.,  dynamically  unstable  configurations)  or  to 
another  apolication  altogether,  the  logic 
associated  with  selecting  the  knot  locations 
should  be  revisited. 


The  second  case  utilizing  straight  line 
segments  splined  together  possessed  only  two 
segments,  n»2  in  Equation  (4),  but  the  knot 
(  5[)  was  allowed  to  be  a  free  variable  and 
determined  by  the  data  reduction  procedure. 
Hence  the  five  unknown  coefficients  become 
CD0>C1 *c2» 61  *nd  cDy  ^see  pi8ure  4). 

Another  set  of  splines  used  in  the  present 
investigation  involves  quadratic  segments 
similar  to  Equation  (2).  Here,  the  knot 
locations  are  free  variables  and  are  determined 
by  the  data  reduction  routine  (see  Figures  5a 
and  5b).  This  expansion  can  be  defined  such 
that  not  only  is  the  Cp  function  continuous  at 
the  knot,  but  also  the  derivative  can  be  made 
continuous  at  the  knot.  Both  of  these 
approaches  sre  discussed.  First,  consider  the 
approach  where  only  Cp  is  continuous  at  the 
floating  knot  location.  This  extension  is 
written  ss 


•••.It  1  ;  1-,  -  -  l.-r.il  *.r  r.vl  .-hr 


I  •  1 
•  • 

I 

I 

I 

I 

_l _ 


Here  CpQ,  Cj,  C2,  '1,  and  Cpv  are 

the  five  free  variables  and  the  function  has  a 

discontinuous  slope  at  Sj  and  S?  are 

the  determining  switches  for  the  polynomials  and 

are  set  similarly  to  those  in  Equation  (4) 

(i.e.,  if  <■  ;^2  then  Sj  *  1  and  Sj  “ 

0  or  if  >‘]2,  then  Si  *  0  and  S-,  = 

This  approach  reauires  the  evaluation  of 
onlv  one  set  of  partial  derivatives  with  respect 
to  C]  or  Cn  depending  on  the  magnitude  of  i’. 

Equation  (5)  can  be  modified  bv  adding  an 
additional  term  such  that  the  slope  of  the  Cp 
vs  curve  is  also  continuous  at  the  floating  knot 
location  to  yield 

v\,  *  c.  *  >:.(.•>  *  V,'?  ‘  * 

„  *  t  1  •  • 


The  unknown  coefficipnts  f  **r  t :  *  express -"in  *r(t 
the  sane  as  for  Equation  (5),  and  the  switches 
are  sec  sinilarly. 
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There  are  manv  other  apline  expansions  that 
could  be  considered  using  five  free  unknown 
coeff icients;  however,  it  is  believed  the  ones 
defined  herein  are  sufficient  to  illustrate  the 
applicability  and  usefulness  of  these 

techniques. 


( 


i  i.  .  i  "v.>  straight  lint  segments  with 
fU’jtir..-  knot  h-i  nr  ion 


b.  CD  vs  8 


5  Two  quadratic  segments  with  floating 
knot  location 


Parameter  Ident if icat ic 


tine  and  distance  data  with  the  numerical 
solution  of  Equation  (1).  This  fitting  process 
is  a  least  squares  technique  and  the  angle  of 
attack  history  is  provided  as  an  input.  ~h.- 
method  used  is  the  one  described  bv  Chapman  an  1 
Kirk?.  This  method  will  be  described  br!*f! 
to  illustrate  the  technique  when  spline  functions 
are  employed  for  expansion  of  the  total  drag 
coefficient.  The  steps  utilized  in  anplvin?  the 
technique  are  as  follows: 

1.  Ident i fv/f ormulate  the  associate! 
equation  of  motion.  For  the  application 
discussed  this  has  been  accomplished  and  is 
Equation  (1). 

2.  Select  an  appropriate  expansion  of  the 
aerodynamic  coefficients.  This  step  js  tin- 
subject  of  this  paper:  the  investigation  o! 
several  expansion  techniques  as  defined  in 
Equations  (2)  through  (61. 

3.  Partially  differentiating  the  equation 
of  motion,  Equation  (11,  with  respect  to  each  of 
the  free  unknown  coefficients,  form  a  set  oc 
parametric  differential  equations.  This  is 
illustrated  below  hv  using  the  expansion  shown  in 
Equation  (51.  Let 


Applying  this  to  Equations  (11  and  (21  the 
following  set  of  parametric  differentia! 
equations  are  derived. 


where  K  is  gA/2m  and  (  )0  implies  the  quantity 
was  evaluated  with  given  coefficients  (either 
initial  guesses  or  corrected  valuesl  at  the  start 
of  each  iteration  cycle. 

A.  Numerically  integrate  the  equation  of 
motion,  Equation  (1),  utilizing  initial  guesses 
for  the  unknown  aerodynamic  coefficients  and 
estimated  initial  conditions  (Xe  and  Xel. 

The  maerical  integration  technique  used  is  a 
Gausa-Nevton  iterative  method. 

5.  Integrate  parametric  equations. 

Equations  (8),  numerically  such  that  the  partial 
derivatives  with  respect  to  each  of  the  free 
unknown  coefficients  are  evaluated. 


The  unknown  free  coefficients  in  the  various  6.  The  method  of  differential  corrections 

^apansions  Equations  (21  through  (6)  are  is  then  used  to  obtain  corrections  to  the  initial 

te rained  by  fitting  the  experimentally  measured  gueasea  of  the  unknown  aerodynamic  coefficients 
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and  estimated  initial  conditions  (Xe  and 
X,,).  This  method  consists  of  expanding  the 
calculated  value  of  position  Xj  caj  about  a 
given  set  of  coefficients  in  a  Taylor  series.  Or 


*  ' : : her  order  Terns 


Note  the  subscript  o  again  indicates  calculated 
values  using  the  given  set  of  coefficients  and 
the  summation,  ignoring  the  higher  order  terms, 
provides  the  contribution  to  Xj  caj  associated 
with  small  changes  in  the  coefficients 
iCj's).  The  sum  of  the  squares  of  the 
residuals,  RSQ,  (difference  between  measured  and 
calculated  downrange  distance  traveled)  is  given 
as 


The  subscript  i  denotes  the  i-th  measurement  and 
i  the  total  number  of  measurements.  Equation 
( 10)  can  accotmnodate  the  reduction/analysis  of 
several  simultaneous  data  sets;  however,  for 
simplicity  che  notation  will  indicate  only  one 
set  of  time  and  distance  measurements  (see 
Reference  A  far  details  on  simultaneous  fitting 
of  multiple  lata  sets  ) . 

Vow  substituting  Equation  ,Q'  into  Equation 
(10),  taking  the  derivative  of  RSQ  with  respect 
to  each  of  the  unknown  coefficients  (Ck)  and 
setting  equal  to  zero,  after  some  manipulation 
one  can  arrive  at  the  following  matrix  equation, 

0  »  f  A  i  k  )  “ 1  Bk  (11) 

where 

i  i  a> 

) 

and 


Here,  C  is  the  matrix  of  corrections  to  be  added 
co  each  of  the  unknown  coefficients  to  be 
letermine!  including  the  initial  conditions  Xe 
and  Xt,.  ',nCp  experimental  data  is  compared 

t>  the  numerically  integrated  position  profile  in 
a  least  squares  sense,  the  corrections  to  the 
unkn'wn  ‘  *  i  c  i  ent  s  and  i-itiil  conditions  can 
’se  lerermineq  using  t_q  tia  r  l  '''1. 

7.  Steps  1  through  *  are  then  repeated, 

isioe  r  he  newlv  ad  iusted  c"e  “  •  c  i  ent  s  and  initial 
conditions,  and  the  process  continues  until 
'osirod  "t'.'.e gence  is  obtained, 

dESVI.TS  AN1>  h' '■''•sRT'oy 


"‘■e  .-sr  mated  total  drag  coefficient  ( Cp ) 

,*■  >ns  obtained  using  various  expansion 

! ecbn 1 ques ,  "quations  (2)  through  (6),  are  shown 
o  r;  cures  s  '"’rough  R.  rach  oc  these  Cp 


expansions  were  evaluated  by  simultaneously 
fitting  four  sets  of  experimentally  measured  time 
and  distance  data  obtained  from  four  separate 
flights  of  a  25m»  spin  stabilized  projectile 
tested  in  the  Aeroballiatics  Research 
Facility*®.  Thia  set  of  data  was  used  to 
evaluate  the  Cp  expansions  because  of  the 
relatively  high  angles  of  attack  experienced 
during  some  of  the  flights  and  the  apparent 
highly  nonlinear  characteristics  of  Cp  with 
total  angle  of  attack.  The  initial  velocity 
(muzzle)  varied  from  3168  to  3245  ft/sec  for  the 
four  flights  and  the  average  mid-range  velocity 
(VREp)  of  all  four  flights  being  3098  ft/sec. 

Here  it  should  be  cautioned  that  because 
some  of  the  expansion  techniques  fit  this 
particular  set  of  data  better  than  others,  it 
doesn't  necessarily  mean  that  one  expansion 
method  is  superior  to  another.  In  fact,  the 
analyst  should  recognize  that  when  selecting  an 
expansion  (whether  it  is  one  of  the  expansions 
discussed  herein  or  another)  for  a  particular 
application  the  inherent  nature  of  the  data 
itself  should  be  the  dominate  consideration.  For 
example,  if  time  and  position  data  obtained  from 
the  flight  of  a  sphere  were  being  analyzed  (angle 
of  attack  is  of  no  concern)  the  expansion  would 
only  include  the  Cp  and  Cpv  terms.  Or, 
if  a  highly  nonlinear  spring-mass-daraper  svstem 
was  being  analyzed,  the  frequency  of  oscillation 
and/or  damping  may  be  modeled  as  a  function  of 
the  displacement  using  multiple  straight  line 
segments  with  or  without  floating  knots.  For 
this  case,  one  mav  possess  information  which 
identifies  where  the  knot  should  be  located. 

Also,  a  good  rule  of  thumb  to  remember  is  to  use 
the  simplest  expansion  which  adequately  matches 
that  particular  set  of  data.  With  these 
considerations  in  mind,  the  Cp  expansions  using 
the  various  techniques  previously  discussed  are 
presented  to  show  the  applicability  and 
versatility  of  these  methods.  Furthermore,  they 
graphically  illustrate  to  anv  potential  user  that 
the  more  conventional  continuous  functions  are 
not  the  onlv  choices  available. 

Figure  6  shows  the  comparison  of  the 
classical  quadratic  (Equation  2)  with  the  fourth 
order  polynomial  including  both  the  ^ a  and 
Cp  terras  (Equation  3).  As  shown  in  this 
figure,  the  fourth  order  polynomial  resulted  in  a 
significantly  better  fit  to  the  four  separate 
flights  than  the  quadratic  (note  sum  of  the 
residuals  squared,  RSQ.  for  both  expansions). 
However,  this  should  not  be  unexpected  since 
the  fourth  order  polynomial  has  three 
additional  unknown  coefficients  (a.  Cp^, 

^f>v),  and  this  added  flexibility  would 
be  expected  to  vield  a  l”wpr  RSO.  Tbe  real  ^ 

question  here  is  whether  or  not  the  dramatic  ri 
in  Cp  at  -he  highest  angles  of  attack  is 
real.  The  four  flights  used  for  these 
comparisons  had  a  total  of  IRQ  data  points  of  ^ 
which  onlv  four  data  points  possessed  angles  o,^ 
attack  greater  than  20  degrees,  '"'once  the  * 
flexibility  mav  he  providing  erroneous  res"  t 
the  higher  angles  of  attack.  This  would  be 
especially  true  if  one  of  the  four  data  point  ^ 
above  2Q  degrees,  is  in  error.  Normally,  beca 
of  the  paucitv  of  data  at  the  higher  angles  o 
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attack,  the  Cp  expansions  would  only  be 
presented  up  to  about  20  degrees.  But  for  the 
purpose  of  this  paper,  the  expansions  are  shown 
0p  to  the  maximum  angle  of  attack. 

ft  is  also  of  interest  to  see  the  effect  of 
the  various  terms  in  the  fourth  order  pol vnoroi a ! 
Equation  CP.  Table  1  presents  the  results 
obtain'*!  fro-’  runn.nc  five  different  c^ses  with 
various  term*-  held  at  zero  in  the  fourth  order 
po l  vnom j  a  1 .  C\i s e  1  shown  in  this  table 
represents  the  classical  quadratic  as  plotted  in 
Figure 


The  comparison  of  the  three  segment  linear 
spline,  Equation  (41  with  n*3,  containing  fix**1 
knot  locations  to  the  two  segment  linear  «p!-r. 
(also  Equation  k  with  n-?)  containing  a  f  !  o  u  i 
knot  location  is  shown  in  Figure  7.  Th : 
illustrates  that  the  fit  containing  three 
segments  is  superior  to  the  two  segment*  (:rr-  5 
with  the  floating  knot  (note  RSQ  for  h'th  tit- 
This  result  is  not  to  he  unexpected  s :  nr**  r!,*- 
qua!  i  tv  of  fit  she*:!  i  inst»rnv«»  hv  incr--j  I  *■  •  »■ 
number  of  segment*.  Ufa  at  is  interest  ’  r 

the  sum  of  the  residual  squared  for  tl-  tSm 
segment  case  is  rl'o  same  as  th**  four:  *»•  ■ 

polynomial  shown  *n  Figure  b  /RSr'  =  .  .  *  ■: 

for  both).  Tt  is  suspected  that  thi*  t  resu 

of  the  fortuitous  knot  location  far'*:  *'-r  ■  . 

bv  dividing  the  f  space  into  three  e.- n  * 
segments)  at  IRP  degrees. 


a 


Case  ■)  represents  the  complete  fourth  order 
polvnomial  expansion,  also  shown  in  Figure  h. 
Cases  2  through  i  represent  the  results  obtained 
when  the  remainina  three  unknown  coefficients, 

C0  ,  Cp^,  and  a,  are  included  in  the 
reduction  routine.  When  viewing  the  results 
shown  in  Table  1,  it  is  apparent  that  the  Cp 
term  is  important  for  this  particular  sec  of  data 
(note  RSQ  for  Case  2  compared  with  Case  1)  and 
appears  to  be  consistentlv  determined  for  all 
cases.  What  is  remarkable  about  the  results 
shown  in  this  table  is  the  dramatic  improvement 
in  the  overall  quality  of  fit  (lower  RSQ)  when 
all  of  the  unknown  coefficients  are  determined 
simultaneously  (see  Case  S).  This  indicates  that 
this  particular  set  of  time-position  data 
requires  a  relatively  complex  expansion  of  Cp 
with  angle  of  attack  and  velocity. 


He.  7  Ch  rr.iri-vT.  iff.  :  t 

si  c-t-r.t  fits  i  st'.  1  ru.,t  : 


a.  Guadratic  splines  plotted  .■ 


H,;.  ‘  (or.parison  of  classic  quadratic  anil 
fourth  order  expansions 


Fig.  8  Cor;  .riser  of  ouadratii  -•  1  *  \  i  f  ■ 

continuous  and  discor.t  in  '.lens 
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The  quadratic  splines  with  a  continuous  Cp 
expansion  only  (Equation  5)  and  a  continuous  Cp 
and  slope  expansion  (Equation  6)  are  compared  in 
Figure  8.  Part  a  of  this  figure  show?  both 
expansions  plotted  against  * ;  whereas.  Part  b 
shows  the  expansions  plotted  vs  1 2 .  Neither  of 
these  expansions  resulted  in  a  fit  to  the  chosen 
set  of  data  as  well  as  the  fourth  order 
polynomial  or  the  three  segment  linear  spline. 

Roih  of  these  expansions  also  contain  a 
condition  in  which  the  coefficient  extraction 
technique  becomes  indeterminate.  This  condition 
exists  when  the  user  permits  the  knot  location  to 
be  a  free  variable  and  the  slope  of  the  first 
segment  is  nearly  equal  to  the  slope  of  the 
second  segment  (Cj-Cj).  When  this  occurs, 
the  minimum  RSQ  is  insensitive  to  the  location  of 
the  knot  and  the  iteration  process  fails  to 
converge.  However,  if  it  is  recognised  that  this 
condition  exists,  the  knot  location  can  be  held 
constant  and  the  iterative  process  again  h*  'omes 
stable.  This  condition  is  also  indicative  of  an 
overlv  defined  svstera.  For  instance,  a  single 
quadratic  expansion  would  fit  the  data  as  well  as 
the  two  quadratic  expansions  splined  together. 

It  seems  for  the  set  of  data  used,  neither 
quadratic  expansion  appears  to  be  advantageous. 
Nevertheless,  this  technique  does  appear  to  hold 
promise  for  modeling  other  aerodynamic 
coefficients  which  are  highly  nonlinear  such  as 
mi  gnus  t*’mi . 

rovri.is  *  N*' 

'irious  total  drag  *oe*f;cienr  expansions 
using  polynomial*  and  '  ir  splines  have  Seen 
ieve loped  and  compared.  A  *»et  >•  four  flights 
was  selected  to  show  the  applicability  and 
versatility  of  these  meth  s. 

AM  of  the  expansions  provided  a  better  fit 
to  the  experimental  data  than  the  classic  linear 
theorv  suggesting  the  need  for  a  relatively 
complex  expansion  of  Cr>  with  angle  of  attack 
and  velocity*  "he  fourth  "r  ler  polynomial  and 
the  three  segment  spline  provided  superior  fits. 
However,  these  results  are  unique  for  this 
'Mrticiinr  lata  set  and  do  not  implv  that  the/ 
are  the  best  expansions  to  use  in  everv  case. 

techniques  as  applied  in  the 
aerodynamic  coefficient  estimation  process  show 
great  potential  for  other  highlv  nonlinear 
anpl i rat  ions.  Moreover,  thev  demonstrate  that 
conventional,  continuous  functions  are  not  the 
'n  \  •/  choices  available  to  the  data  analyst. 


many  of  the  program  modifications  required  during 

this  investigation. 
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